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FIGURE 12 

Pseudoscalar decay constants from Ref. 16 (quenched) and from Ref. 26 (burst), 
at sea quark mass aniq = 0.01 and equivalent lattice spacing. 
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FIGURE 11 

Lattice 1/ fy from Ref. 26. The labeled points are physical particles Results from 
simulations with sea quark mass in lattice units aniq = 0.01 are shown in squares, 
and for sea quark mass aniq = 0.025 in diamonds. 



have to fight your way to the numbers. 
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FIGURE 10 

Argus data for the B ^ C form factor from Ref. 25 in circles, with the lattice 
results of Ref. 22 as crosses. (This figure is Fig. 3 of Ref. 22.) 

Simulations remain unwieldy. Doing almost anything requires at least a year of work on 
a supercomputer. This project length seems to be an invariant-as computers improve, our 
standards have gone up. Nevertheless, the continued improvement in computer hardware 
allows us to tackle more and more complicated (interesting?) projects, so that the field 
will continue to advance even in the absence of new ideas. 

The most interesting new ideas, which might lead to improved calculations on smaller 
computers, are concerned with the question: Can one find a more complicated discretiza- 
tion which allows one to work at bigger lattice spacings? Doubling the number of terms in 
the lattice action roughly doubles the amount of work, while halving the lattice spacing at 
fixed simulation volume increases the work by a factor of 16. This subject is under active 
study?^ 

The best way to end the talk is to quote a previous speaker (Prof. B. Frois): "You 
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FIGURE 9 



Lattice data from various combinations of heavy and light quark masses from Ref. 
22 with some theoretical curves superposed. (This figure is Fig. 1 of Ref. 22.) 



SUMMARY 



Present day lattice calculations are able to produce ten to fifteen per cent numbers for 
a wide variety of physical observables. Most of the uncertainties are systematics limited 
(at the cost of large amounts of computing to beat down statistics). The major systematic 
is the lattice spacing. It is just not understood how small the lattice spacing should be so 
that lattice calculations are insensitive to it (or more precisely, so that all physics on scales 

(\>-7 

less than a are perturbative). There are claims that heavy quark spectroscopy only needs 
a ~ 1/5 fm. However, glueball spectroscopy at that lattice spacing shows a dependence. 
The D-meson decay constant needs a ~ 0.08 fm. The minimum lattice spacing is probably 
process dependent. 
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Table I. Predictions for heavy-light decay constants from Ref. 16, showing various 
uncertainties. 



particle 


fp, MeV 


fitting/extrap 


scale 


large- am 


B 


187(10) 


±12 


±15 


±32 


Bs 


207(9) 


±10 


±22 


±32 


D 


208(9) 


±11 


±12 


±33 


Ds 


230(8) 


±10 


±18 


±28 



line), a fit to the lattice data (solid line), and one theoretical prediction of the curve [24] 
(dotted curve). Fig. 10 shows a comparison with the experimental data from ARGUS^^. 

Of course this is just the beginning of these calculations. The lattice can be used to test 
the universality of the Isgur-Wise function at the B and D-meson masses, and determine 
the dependence of corrections to it on the heavy quark mass. 

Testing the Quenched Approximation 

All of these calculation are performed in the quenched approximation. There is an 
unknown systematic associated with throwing away the sea quarks. The only way I know 
to really test it is to repeat the simulations with dynamical sea quarks. That is very 
expensive. However, there are a few tests already in the market. Fig. 11 shows the decay 
constant of vector mesons parameterized by 

{V\V^\0) = -j^mle^. (20) 

from a simulation with Wilson valence quarks and two flavors of dynamical staggered 
quarks [26]. The two plotting symbols are for two different values of the sea quark mass 
(the lattice spacing is about 1/a ~ 2 GeV). Clearly the effects of sea quarks are small. 
As a second example, we display the pseudoscalar decay constant from the same data set 
against the results of Ref. 16 in Fig. 12. Our data (the burst) is at about the same lattice 
spacing as the quenched data plotted as squares. If there is an effect of sea quarks, it is 
not very large. 
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FIGURE 8 

/7\^a/M/(1 + as/7rlog(ma)) vs. 1/M, from Ref. 16. The results from two lattice 
spacings (/3 values) are shown, to give a rough idea of lattice spacing systematics. 



(Cjg is a short-distance perturbative factor). One would like to calculate i^o(^ '^O from first 
principles. While there has been some discussion of how to do this with infinite mass heavy 
quarks on the lattice^\ another technique is just to calculate a form factor on the lattice 
and fit it to the form of Eqn. (18). Bernard, Shen, and Soni^^ have recently published a 
preprint which does just that, by measuring the form factor 



{D{v')\cjiyc\D{v)) = mj:,Ccc{i^)Co{v' ■ v, ij){v + 



(19) 



The lattice calculation spans the range Q < v ■ v' < 1.2 while real- world data ranges over 
1.1 < V ■ v' < 1.5. 

Their results are shown in Figs. 9 and 10. Fig. 9 shows lattice i^o(^ '^O ^ function of 
{v -v' ) for a variety of light and heavy quark masses (see Ref. 22 for details) and also shows 
various theoretical predictions: an upper bound on the Isgur-Wise function [23] (dashed 
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A Case Study-Heavy Meson Decay Constants 



The decay constant of a pseudoscalar meson M is defined as 

(0|V^7o75V'|M) = fu^M- (16) 

Decay constants are interesting because some of them (vr and A') are measured and provide 
a benchmark for lattice calculations, while some of them are not measured and allow 
predictions (D, Dg, and B). They probe very simple properties of the wave function: in 
the nonrelativistic quark model 

/„ = 4a (17, 

where tp{0) is the qq wave function at the origin. For a heavy quark (Q) light quark (q) 
system tp{0) should become independent of the heavy quark's mass as the Q mass goes to 
infinity, and in that limit one can show in QCD that ^/rn]^fJ\^ approaches a constant. It 
is believed that CP nonconserving amplitudes are proportional to and so knowledge 
of provides information about CP nonconservation in the B system. 

One way to compute the decay constant is to put a light quark and a heavy quark on 
the lattice and let them propagate. It is difficult to calculate directly on present day 
lattices because the lattice spacing is much greater than the b quark's Compton wavelength 
(or the UV cutoff is below m^). In this limit the b quark is strongly affected by lattice 
artifacts as it propagates. However, one can make infinite on the lattice and determine 
the combination ^/rf^gfJ^ in the limit. Then one can extrapolate down to the B mass and 
see if the two extrapolations up and down give the same result. Until a year or so ago the 
two methods did not give consistent numbers. However, the present situation is that one 
can reliably compute fjj and in quenched approximation. 

As an example, results from a recent calculation by Bernard, Labrenz, and Soni^^ is 
shown in Fig. 8. What is plotted is fj^\/M /{I + ag/vr log(Ma)) vs. 1/M; the extra term 
is a perturbative correction to the static heavy quark formula. Removing it allows one 
to interpolate to infinite quark mass. Their predictions are reproduced in Table I. There 
are several other lattice predictions of these numbers. (See Ref. 17 for a compilation.) 
They differ in detail, but all give numbers in the range of those of Table I. There are two 
experimental measurements of Z^,^. They are 232±45±20±48 MeV [18] or 344±37±52±42 
MeV[19]. The error bars are too big for a serious comparison. 



The Isgur-Wise Function 



The physics of systems containing a heavy quark and a light quark has a very simple 
limit as the mass of the heavy quark goes to infinity. The physics of the light quark becomes 
independent of the mass or other properties of the heavy quark. (For an extensive review, 
see Ref. 20.) In particular, the form factor in B ^ D semileptonic decay is described 
by a universal function called the Isgur-Wise function i^o(^' ' ^7 A*) which depends on the 
four- velocities of the two heavy quarks: 

{D{v')\c-ij,b\B{v)) = V"^5"^L>C'c&(/^)Co(^^' • V, iJ.)(v + v')i, (18) 
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IV. MATRIX ELEMENTS 



Most of the matrix elements measured on the lattice are expectation values of local op- 
erators composed of quark and gluon fields. The mechanical part of the lattice calculation 
begins by writing down some Green's function which contains the local operator (call it 
J{x)) and somehow extracting the matrix element. For example, if one wanted (0| J(a;)|/i) 
one could look at the two-point function 

Cjo(t) = 5](0|J(x,t)O(0,0)|0) (10). 

X 

Inserting a complete set of correctly normalized momentum eigenstates 



and using translational invariance and going to large t gives 



1 ^ 



Cjoit) = e — . (12) 



A second calculation of 



Cooit) = 5](0|O(x,t)O(0,0)|0) = (13) 

X ^ 

can be used to extract (0| J|y4.) (fit two correlators with three parameters). 
Similarly, a matrix element {h\J\h') can be gotten from 

CABit.t') = J](0|O^(t)J(x,t')O5(0)|0). (14) 

X 

(Can you see how?) 

These lattice matrix elements are not yet the continuum matrix elements. The lattice 
is a UV regulator and changing from the lattice cutoff to a continuum regulator (like MS) 
introduces a shift 

(/|0-"^(/. = l/a)\z)j^ = a^il + g(C^- Q„«) + • • .)(/|0'««(a)|z) + 0(a) + .... (15) 

The factor converts the dimensionless lattice number to its continuum result. The 

0(a) corrections arise because the lattice operator might not be the continuum operator: 

df I dx = {f{x -\- a) — f{x))l a -\- 0{a). The C's are calculable in perturbation theory. There 

are a number of tricks/deep theoretical ideas for achieving a more convergent perturbation 
1 

expansion. 
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FIGURE 7 

Glueball masses from the UKQCD collaboration [14] The black symbols are 
measurements and the open symbols are just lower limits on the masses. (This 
figure is Fig. 1 of Ref. 14.) 

Glueballs 

People have been calculating glueball masses in QCD for many years. Unfortunately, 
all the high statistics calculations use quenched approximation; you would find nothing 
interesting in any of the calculations which include sea quarks. A representative recent 
compilation"'^'^ is shown in Fig. 7. All of these states are above vrvr threshold and so one 
would worry how these numbers will change when sea quarks are included. That's an open 
problem. 
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FIGURE 6 

The strong coupling constant at the mass of the Z. The lattice number is shown 
along with other determinations from various experiments, from Ref. 13. 



the other ones. While I feel that the lattice number will eventually turn out to be more 
reliable, since it only uses spectroscopy as an input, no jet physics, I don't think that is 
the case yet. I would really prefer to see checks at smaller lattice spacing and (eventually) 
the full lattice simulation done in the presence of dynamical fermions. Maybe I just don't 
know where all the bodies are buried in the e~^e~ analyses! 

Heavy quark systems are better than light quark systems for exploring the effects of 
quenching, because light quarks only modify the potential between the heavy quarks. One 
can play with models to understand their effects. In contrast, we don't really understand 
the effects of sea quarks on light hadron spectroscopy. For example, why are the rho and 
omega mesons nearly degenerate even though the decay width of the rho is so much greater 
than the omega's? 
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FIGURE 5 

Masses of the four P-wave cc states, from Ref. 10. 



Next, one must determine the couphng constant. This is done through the short range 
QQ potential, 

n?') = ^«.(?')^ (8) 

where g ~ vr/a. Lattice perturbation theory does a good job at the lattice spacings of the 
simulation, and so one can measure the short distance potential, extract on the lattice, 

and carry out the conversion. Doing so gives (^j^ (5GeV) = 0.140(4). 

The main problem with the calculation as it presently stands is that it is done in 

quenched approximation. One must somehow convert a— ^ to a— ^ . The authors of 
^ _ _ MS MS 

Ref. [11] do this by running the coupling constant down from the upsilon mass to the typical 
bound state Q with no flavors, then out with four flavors. They flnd that this shifts 
by 25 lb 6 per cent-that is, the uncertainty in the amount of the shift is itself 25 per cent. 

71 f — 4 

This seems conservative. The bottom line is that a— — (5 GeV) = 0.174(12). A second, 

n f — 4 

completely separate calculation using nonrelativistic quarks [12] gives (^-^^g (5 GeV) = 



0.170(10). Finally, running down to the Z-mass gives 



a 



nj=A 
^MS 



(Mz) = 0.105(4) (9) 



This result is shown with other determinations of ag in Fig. 7. It is a little bit lower than 
most of the results from LEP. 

At the conference, several theorists wanted to make a big stir that the lattice number 
is somehow better than the other determinations, and that it should be trusted more than 
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Edinburgh plot from simulations with two flavors of dynamical staggered fermions. 
Data are by Bernard, et. al.[6] (cross and fancy square), Bitar, et. al. [7] (square 
and diamond), F. Butler, et. al. [8] (fancy diamond) and M. Fukugita, et. al. 
[9](fancy cross and burst), while the upper curve is a theoretical prediction in 
strong coupling (from Ref. 5) and the lower curve is the extrapolation of the data 
at /3 = 5.7 to zero quark mass. 



nect something perturbative to lattice perturbation theory and then to a continuum (MS) 
number. The calculation has been done by two groups: El-Khadra, Hockney, Kronfeld, 

11 • 1 

and Mackenzie , and Davies, Lepage and Thacker . 

The calculation begins by noticing that the mass difference between the lightest S-wave 
and P-wave QQ mesons is nearly independent of quark mass (in the tjj system it is 460 
MeV, in the T, 430 MeV. Since this difference is independent of the quark mass, one does 
not have to tune the quark mass on the lattice in order to measure it. So the lattice S-P 
mass splitting gives the lattice spacing. 
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FIGURE 3 

Ratio of lattice masses to the rho mass, after extrapolations to infinite simulation 
volume and zero lattice spacing, from Ref. 4. Circles are real world data, squares 
from simulations 



Heavy Quarks 



People have been using the lattice to do calculations of heavy quark systems, too. As 
an example. Fig. 5 is a picture of the fine structure splitting in charmonium, by me and 
M. Hecht"'^'^. It looks just like the wallet card, although with large errors. 

The most interesting application of heavy quark physics recently reported is the cal- 

The idea is that the lattice needs one parameter to set the 



culation of Ai 



or a-, 



MS MS- 

scale. The easiest parameter to determine on the lattice is a mass or mass difference and 
the idea is to use a mass difference to find ^j^- The calculation has two parts. First, 
one must do a long distance calculation to find a mass; the lattice spacing comes from 
a = Mffa/Mff(expt). Next, one must do a short distance calculation on the lattice to con- 
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Edinburgh plot prepared by me from the data of Ref. 4, showing ratios at several 
values of the lattice spacing (different /3's). The octagon shows the expected result 
at infinite quark mass, and the question mark is the real world value. 



As an example, I show In Fig. 4 an Edinburgh plot for simulations with two flavors of 
dynamical staggered fermions. I have connected the points with the same lattice spacing 
(same fi). Again, the iV/p ratio appears to "settle" a bit as fi increases. The upper curve 
is an analytic calculation at inflnitely strong coupling, /3 = oo, where the lattice spacing 
is about 1/2 Fermi [5]. The lower curve is the extrapolation of the fi = 5.7 data to zero 
quark mass. We see that in that limit the iV/p ratio is still too large. The lattice spacings 
here are all much larger than in Fig. 2-at fi = 5.445 a = 0.22 fm, at /3 = 5.6 a = 0.11 fm, 
and at /3 = 5.7 a = 0.089 fm. (The relation between fi and the lattice spacing is different 
at = then at 7^ because the QCD /3-function depends on nj.) These calculations 
have a ways to go. 
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FIGURE 1 

(a) A typical correlator showing good exponential fallofF (the correlator has peri- 
odic boundary conditions in the time direction), (b) Feynman diagrams for meson 
and baryon correlators. 



about 0.14 fm) to the smaller lattice spacing (/3 = 5.93, 6.17, a down to about 0.07 fm). 
The authors of Ref. [4] have extrapolated their masses in a and L and present the limits 
in Fig. 3, as a plot of mass divided by Mp at M^j- = 0. The agreement with observation is 
spectacular. 

All of these calculations are done in quenched approximation. Simulations with dy- 
namical fermions are much more expensive and the data is correspondingly more meagre. 
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III. SPECTROSCOPY 



All lattice calculations begin with spectroscopy. In order to measure the mass of a 
hadron which has some set of quantum numbers, invent an operator J which has the same 
set of quantum numbers and compute 

C(t) = (0|J(t)J(0)|0). (4) 

Using the Euclidean version of the Heisenberg equation of motion 



and inserting a complete set of energy eigenstates, we find 



which at big t goes over to 



J{t) = exp{Ht)J exp{-Ht) 

C{t) = J]|(0|J|n)|2exp(-E„t) 

n 

C{t) ~ |(0|J|l)|2exp(-Eit) 



(5) 
(6) 

(7) 



where Ei is the lightest state with the quantum numbers of J. The exponential falloff of 
the correlator gives us the mass, while its intercept gives us a matrix element (0|J|1). 

For bound states of quarks the operator C{t) is basically the Feynman graph shown in 
Fig. 1: it is made of the appropriate number of quark and antiquark quarks propagating 
in the background of gluon fields in your album of configurations. 

In the old days (pre-1988) the operators J were local currents, like V'TsV' for the pion. 
Nowadays we use some big extended operator like cf){x ,y)il){x)^^il){y) which "looks 

like" a hadronic wave function. Then the computer has to do less work to filter out the 
lightest state. This means that Eqn. (6) takes its asymptotic form Eqn. (7) at a small t 
value, while the signal is still large. 

Spectroscopic studies in QCD involve light quarks, heavy quarks, and glueballs, so I 
will say a few words about each. 



Light Quarks 



Generally in lattice calculations people try to deal with dimensionless quantities as 
much as possible, since they are independent of the precise value of the lattice spacing. 
In spectroscopy, people present their data on so-called "Edinburgh plots," M]\[/Mp vs. 
Mt,/Mp or "APE plot," Mj^/Mp vs. {MT^/Mpf. (The names are after the collaborations 
which invented the plots). 

The most interesting recent quenched calculation is by a group from IBM which built its 
own computer to do QCD [4]. Fig. 2 shows their data plotted by me on an Edinburgh plot. 
There appears to be a small change between the data at larger lattice spacing (/3 = 5.7, 
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Brillouin zone pa = vr. With staggered fermions these extra states are treated as extra 
spin or flavor degrees of freedom. 

5. "Wilson fermions" add extra terms to the action to raise the energy of the p = ir/a 
modes and ehminate the extra degeneracy. The physical quark mass must be derived 
from the simulation in terms of an input parameter called the hopping parameter k, 
and a measured parameter Kc'. aniq = 1/2(1/ k, — l/nc). Wilson fermion calculators 
like to use k, in their graphs. 
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typically performs a calculation at an unphysical value of the light quark mass and 
then tries to extrapolate to niq = 0. 

5. Sea quarks are a problem because of Fermi statistics, which effectively introduces 
long range interactions among the quarks. There are techniques for dealing with this 
problem^'*^ but they make QCD with dynamical fermions orders of magnitude more 
difficult than if the sea quarks were not there (and the difficulty scales inversely as a 
power of the quark mass). A rather drastic approximation called the quenched approx- 
imation neglects this problem simply by throwing away all the sea quarks. This is an 
uncontrolled approximation which people do mainly because the alternative (keeping 
light sea quarks) is too time consuming for the computer. 

All these constraints add up to a very hard numerical problem. We use the fastest 
supercomputers available. Cray's are usually too slow. Some groups have built their own 
computers. One of the projects I belong to used half of a Connection Machine CM-2 (at 
a speed of about 3 1/2 Gflops) for about two years. This is not considered an excessive 
amount of resources. 

Finally, there are two more problems to watch out for. 

6. In a lattice calculation all observables are measured on the same set of lattices and are 
highly correlated. There are methods for dealing with correlated data. Some lattice 
practitioners use them. My advice is that if the paper you are reading does not make 
some attempt to deal with the correlations which are present in its data (or is not 
aware that its data is correlated), you should discard the paper. 

7. The major problem facing lattice calculations these days are systematics: Quenching, is 
a small enough, is L big enough, is the quark mass small enough? Lattice calculations 
produce as output not a hadron mass niff but the combination amjj. One finds a by 
dividing amjj by a measured mjj (in MeV). The problem is, which mass to use? Most 
lattice calculations only reproduce mass ratios at the ten or fifteen per cent level, so the 
lattice spacing is uncertain at that level. This uncertainty propagates into essentially 
all interesting calculations. 

I would be remiss if I did not provide you with a small glossary of lattice terms in 
order to enable you to read the literature: 

1. Lattice people define /3 = 6/(/^ where = (/^/47r. Here g is the color coupling constant 
measured at a momentum scale Q ~ vr/a, so bigger fi corresponds to smaller a. 

2. "Link"- The vector potential A^{x) has an orientation and so instead of being defined 
on the sites of the lattice, is defined on the links joining adjacent points x and x + an^. 
For technical reasons lattice people use the "link variable" U^{x) = ex.]){igaA^{x)) in 
simulations rather than the vector potential. 

3. "Plaquette"- The lattice analog of the gauge action F^j^(x) is the product of four links 
about a unit square or "plaquette" on the lattice. 

4. "Staggered fermions"- On the lattice the quark energy momentum dispersion relation 
changes from its continuum value E = p + m to sinh Ea = sin pa + ma. This 
has low energy states near pa = and degenerate extra solutions at the ends of the 
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II. HOW LATTICE CALCULATIONS ARE CARRIED OUT 



Lattice calculations are performed using the Euclidean path integral formulation of 
quantum field theory. If we have some field theory with field variables (f> {(f> could be 
quarks, gluons,. . . ) and a Lagrange density C{cf)), we define an analog of the partition 
function in statistical mechanics 

(here = {x^it)). The expectation value of any observable 0{cf)) is given by 



(O) = [#(x,t)]0(<^)exp(- j d''xC{<P)). (2) 

To be able to perform calculations in any quantum field theory one must introduce a 
short distance cutoff which regulates the ultraviolet divergences. We do that by replacing 
continuous space time by a lattice of grid points x = ax{ where a is the lattice spacing, and 
defining the field on those grid points (f>{x) — )• (f>{ = (f>{x{). Then the functional integrals 
Eqns. (1) and (2) become ordinary integrals of very high dimensionality. One evaluates 
Eqn. (2) using importance sampling: somehow one creates an album of snapshots of the 
field variables (f>{ where the probability that a particular configuration is present in the 
album is P{<f>j) = exp(— ^i^d then 

1 ^ 1 

(^) = ^E^(^.) + o(:^). (3) 



The generation of the album is done using Monte Carlo techniques not too different in 
principle from the ones you would use in an experiment to generate Monte Carlo events. 

Lattice calculations are hard for several reasons: 

1. The lattice spacing should be small-small enough that physics on a size scale less than 
a lattice spacing can be described using perturbation theory. 

2. The size of the simulation volume L'^ should be greater than the physical size of the 
hadrons. This point is in conflict with item 1. The number of grid points is n = (L/a)^. 
A gluon fleld is a three by three complex matrix per each direction on each lattice site, 
or 72 real numbers per lattice site. Fermions have four spins and three colors or 24 real 
numbers per site. Typical simulations have lattice spacings around 1/10 fermi (within 
a factor of two) and a number of mesh points ranging from 16^ X 32 to 24^ x 40 to 32"*: 
the end is not yet in sight! 

3. One needs a lot of statistics-tens to hundreds of uncorrelated lattice measurements. 

4. It is very hard to compute with light (u,d) quark masses at their physical values. 
On the lattice calculating a quark propagator Gq{x,x') involves inverting the matrix 
problem (Ip — m)Gq(x, x') = S^(x — x'). The matrix becomes singular as niq — )• 0. One 
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I. INTRODUCTION 



This lecture is an introduction to lattice calculations in quantum chromodynamics 
for the non-expert "consumer." Lattice methods are presently the only way to perform 
calculations of masses and some matrix elements in the strong interactions beginning with 
the Lagrangian of QCD and including no additional parameters. By "consumer" I mean 
a person who might want to use a lattice calculation (as an input to a phenomenological 
calculation or to compare to her experiment, for example), and is not really interested in 
doing the lattice calculation herself, but would like to be able to judge the reliability of 
calculations in the literature. 

There are many good reviews and introductions to lattice gauge theory and its use 
in QCD?^ The lattice community has a large annual meeting and the proceedings of those 
meetings (Lattice 'XX, published so far by North Holland) are the best places to find the 
most recent results. However, as in any large community with its own set of problems, 
most of the papers in those proceedings tend to talk to each other in a language which 
is rather opaque to nonmembers. My goal is an impressionistic overview of the field as it 
presently exists, which might be useful to an outsider. 

The bottom line is that for the past one or two years there have been a lot of lattice 
calculations of masses and matrix elements which agree with experiment at the ten to 
fifteen per cent level. 

I will begin with a very superficial overview of how lattice calculations are performed. 
Then I will turn to a set of case studies: spectroscopy of light hadrons, of heavy quark 
systems, and of glueballs, then two case studies of matrix elements: the decay constants 
of D- and B-mesons, a recent calculation of the Isgur-Wise function, and some pictures of 
the effects of sea quarks on simple matrix elements. 
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